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Abstract 

The phase structure of the generalized Yang-Mills theories is studied, 
and it is shown that almost always, it is of the third order. As a specific 
example, it is shown that all of the models with the interaction y\ (rij — 
j + N) 2k exhibit third order phase transition, (nj is the length of the j-th 
row of the Yang tableau corresponding to U(JV).) The special cases where 
the transition is not of the third order are also considered and, as a specific 
example, it is shown that the model ^T( n j — J+-^0 2 + 5 ^2j( n 3 ~ j + N) 4 
exhibits a third order phase transition, except for g — 27 'it 2 / '256, where 
the order of the transition is 5/2. 
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1 Introduction 



Recent works by many authors suggest that the pure Yang-Mills theory in two 
dimensions is equivalent to a string theory [1-7]. The reason to investigate 
such a simple theory in such great detail is partly to shed light on the strong- 
coupling limit of pure QCD in four dimensions. The string picture of the two- 
dimensional Yang-Mills (YM2) is also interesting on its own, as an example for 
nonperturbative analysis of a quantum field theory. 

The starting point for making the correspondence of YM2 with string theory 
is the study the large- N limit of YM2. For example, as it is shown in jjj, ||, and 
[Q, a gauge theory based on SU(iV) is split at large N into two copies of a chiral 
theory, which encapsulate the geometry of the string maps. The chiral theory 
associated to the Yang-Mills theory on a two-manifold M. is a summation over 
maps from the two-dimensional world sheet (of arbitrary genus) to the manifold 
A4. This leads to a 1/N expansion for the partition function and observables 
that is convergent for all of the values of area x coupling constant on the target 
space A4, if the genus is one or greater. 

It is well known that YM2 is defined by the Lagrangian tr(F 2 /4) on a Rie- 
mann surface, where F is the field-strength tensor. This theory is equivalent 
to the so-called B—F theory, characterized by the Lagrangian itr(BF)+tr(B 2 ), 
integrated over the auxiliary field B. This theory is essentially characterized 
by two important properties: the lack of propagating degrees of freedom, and 
invariance under area-preserving diffeomorphisms. This theory, however, is not 
unique in possessing these two characteristics. These properties are also shared 
by a large class of theories, called the generalized two-dimensional Yang-Mills 
(gYM2) theories. These theories are defined by replacing the term tr(B 2 ) in the 
YM2 Lagrangian by an arbitrary class function of B || . 

Many of the physical features of gYM2 have been studied. Several aspects 
of the partition function of these theories has been discussed in [8-10]. The 
generating functional of the field strength for arbitrary two-dimensional ori- 
entable and non-orientable surfaces has been calculated in |ll[ (and in [jl2[ |l3| 
for the YM2's). In all of these calculations, the solutions appear as some infi- 
nite summations over the irreducible representations of the gauge group U(iV) 
(or SU(iV)). In the large-TV limit, however, these summations are replaced 
by suitable path integrals over continuous parameters characterizing the Young 
tableaux, and a saddle-point analysis shows that the only significant represen- 
tation is the classical one, which maximizes some effective action. 

The continuous function characterizing the representation is a constrained 
function, as the length of the rows of the Yang tableau is nonincreasing. This 
means that the problem of minimizing the effective action is a constrained one, 
and in fact the constraint is a nonholonomic one. For small values of the surface 
area, the classical solution satisfies the constraint; for large values of the sur- 
face area, it does not, and the dominating representation is not the one which 
minimizes the effective action. This introduces a phase transition between these 
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two regions, which has been studied in (for example) [jli], [f5] for YM2. This 
same problem has been studied for some special cases of gYM2 's Q . In |l7j , 
the density function p has been obtained for an arbitrary gYM2 in the weak- 
coupling region, p is essentially the density function of the length of the rows of 
the Young tableau, and by weak-field region it is meant that the surface area A 
is less than some critical area A c . Moreover, calculating p in the strong region 
A > A c (near the critical point A c ), it is shown that the phase transition of 
the theory corresponding to the interaction J2j( n j ~ 3 + ^0 4 i s °f the third 
order, the same behavior seen for YM2. Finally in |{is| , the phase transition of 

Ej ( n 3 -J+ N ) 6 and Ej ( n 3 ~3+ N Y + 9( n i ~3 + ( for 9 « !) has been 
proved to be of the third order. In all of these works, the procedure of investi- 
gating the phase transition is based on the direct calculation of p s (the density 
function in the strong region A > A c ). The form of this function is complicated 
and obtaining a closed form for it, for arbitrary gYM^'s, is not easy, even near 
the critical point. 

In this paper we are going to investigate the order of the phase transition 
for a typical gYM2, by a different method. In this method, only the behavior 
of p w (the density function in the weak-coupling region) is needed, and in fact 
this behavior is needed only near the critical point. In this way, one finds 
the phase transition of almost all of gYM^'s is of the third order. Nontypical 
gYM2's are also investigated and the criteria for a theory to be typical, that is 
its phase transition be of the third order are also studied. It is shown that if the 
second derivative of the weak-region density function is negative at its absolute 
maximum for A = A c , and if its first derivative with respect to A is positive 
at the same point, then the theory exhibits third order phase transition. It is 
further shown that theories with monomial interactions satisfy these criteria. It 
is also shown that if the second criterion is violated, then the order of transition 
will be an integer times 3. Finally, a theory with quadratic as well as quartic 
interaction is analyzed and it is shown that it is typical except for a specific 
value of the coupling constant, where it exhibits a transition of the order 5/2. 

Let us investigate this nontypical behavior more closely. It was stated that 
if, at the critical point, the first derivative of the absolute maximum of the 
density function with respect to A is positive, and the second derivative of the 
density function at its absolute maximum is negative, then the theory is typical, 
that is the order of the transition is 3. What is known, is that at the critical 
point the absolute maximum of the density function is an increasing function of 
A. The reason is that for A < A c , this absolute maximum is less than one, as 
the system is in the weak region. For A > A c , this absolute maximum should 
be greater than one. So it should pass one at A = A c . This shows that the 
first derivative of this absolute maximum is nonnegative, but it may be zero at 
A = A c . Theories for which this derivative is zero exactly at A — A c , are fine 
tuned. Normally, this derivative is positive. The second criterion for a theory 
to be typical can also be discussed the same way. At the absolute maximum 
of the density function, the first derivative of the density function is, of course, 
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zero. The second derivative should be nonpositivc. Theories for which this 
second derivative is zero exactly at the transition point, are fine tuned. Other 
theories are typical and exhibit a third order phase transition. In fact, in the 
example which shows a transition of the order (i.e. exponent) 5/2, exactly this 
criterion is violated. The parameters of a theory are tuned so that at the critical 
point, the second derivative of the density function is zero, and the first nonzero 
derivative of that function is the fourth derivative. 

The scheme of the paper is the following. In section 2, after a review of 
gYM^'s, the criteria for the existence of the third order phase transition is 
studied and the order of the transition for nontypical theories is also studied. 
In section 3, it is shown that theories with monomial interaction are typical. In 
section 4, a theory with a quadratic and quartic interaction is studied, and it is 
shown that it exhibits third order phase transition except for a specific value of 
the coupling constant. It is shown that for this value, the order of the transition 
is 5/2. Finally, in section 5 the standard technique is used to investigate the 
critical behavior of this theory, and again it is checked that the order of the 
transition is 5/2 at this specific value of the coupling constant. 

2 Typical gYM2 theories 

The partition function of the gYM2 on a sphere is [jlo| [nj 

r 

where r's label the irreducible representations of the gauge group, d r is the 
dimension of the representation r, A is the area of the sphere, and A(r) is : 

A M = £^T^(r), (2) 

k=l 

in which C k is the fc'th Casimir of group, and afc's are arbitrary constants. 
Now consider the gauge group U(iV) and parameterize its representation by the 
integers n\ > n-i > ■ ■ • > n^- It is found that 

* - n 

l<i<j<N v ■' 7 

N 

C k = Y}(^+N~i) k ~(N-i) k ]. (3) 

For the partition function (1) to be convergent, it is necessary that p in (2) be 
even and a p be positive. 
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In the large- N limit, the above summation is replaced by a path integration 
over the continuous function 

<j>{x) := -n(x) - 1 + x, (4) 

where 

0<x:=i/N <l and n(x) :— rn/N. (5) 
The partition function (Q) then becomes 

Z= ( J] d<j)(x)e s ^\ (6) 
J a<x<i 

where 



5(0) = i-A J dxG[<f>(x)} + J dx J dy log \<j>{x) - , (7) 

apart from an unimportant constant, and 

G(0) = £(-1)W- (8) 
fc=i 

Introducing the density function 

*<*>! = s&y < 9 > 

it is seen that it satisfies 

/ p(z)dz = l, (10) 

Jb 

where [b, a] is the interval corresponding to the values of 4>. If G(cj>), and therefore 
p(z), is even, then b = —a. The condition n± > > ■ ■ ■ > demands 

p(z)<l. (11) 

As N — > co, only that representation contributes to the partition function (Jfy 
which maximizes S. This representation is found to be specified by the solution 
of the following integral equation M . 

g{z) = p rpjz^ t (12) 
Jb z - z 

where P indicates the principal value of the integral, and 

g(z) := ^G'(z). (13) 
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The free energy of the theory is defined through 



from which we have 

F'(A) = f dz G{z)p{z). (15) 

Jb 

Using the standard method of solving the integral equation (|lj), the density 
function p(z) is obtained in terms of the two parameters a and 6, where o and b 
themselves must be calculated by solving two other equations |l7| . If G is even, 
the results are 



and 



= V^ZZ y ^~ 1)u » 2 "zV 2n+g+1) (o), (is) 

I K ' 7T ^ 2 n n\(2n + q + 1)! y w ' K ' 

n,q=0 y 



(2n 1)" an (an -i) (Q) = x 
2 n n!(2n- 1)! y w v ' 



n=0 



where is the nth derivative of <?(#). As it is seen from (p7[), one cannot find 
the closed form of a for arbitrary gYM 2 's (that is, for arbitrary G's). 

The function p(z) found from (|l6|), obviously depends on the area A. Its 
absolute maximum is an increasing function of A. But the constraint (^Tj) tells 
that this maximum must not exceed 1. So, the function p(z) obtained from ( |l6| ) 
is valid only for A less than some critical value A c . A c is the value of A at which 
the maximum of p becomes 1. The region A < A c is called the weak-coupling 
region. In the region A > A c (the strong-coupling region), the following ansatz 
for p(z) is used ||l7f . 

*<«>-{«.>, (18 » 

where L' is some interval containing the point (or points) maximizing p. The 
form of p s is usually found to consist of some integrals not expressible in terms of 
elementary functions. In special cases, it has been calculated perturbatively for 
areas slightly more than the critical area for example). Having 

obtained p s , the difference between the free energies in strong and weak regimes 
(F s — F w ), for areas slightly more than the critical area is calculated to obtain 
the order of the transition. 

Consider this same problem of phase transition from a different point of 
view. Denote by p w {z) the weak-region density function, and by p s {z) the 
density function (|T§|). Suppose that p w (z) attains its (absolute) maximum at 
z = zq. We have 

Pw{z) = Pw(z ) + ^p'u(z )(z - Zq) 2 H , (19) 
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in which p^(zo) is assumed to be nonzero. For A slightly more than A c , p w (zo) 
is slightly more than 1. We want to find the points z\ satisfying p w (zi) = 1. It 
is easy to see that 



\zi - Zo\ « y/2[l- p w (z )]/p^(z Q ) ~ a 1/2 , 
where a := p w (z<j) — 1. Next, consider F s — F w : 

6F 



(20) 



F[ Ps ] - F[ Pw ] 



Sp(z) 



[Ps{z) - p w {z)]dz 



6 2 F 



6p(z)5p(z<) 



[p s (z) - p w (z)} 



x[p s (z') - p w (z')]dz dz + 
As p w is the classical solution, it satisfies 

s 4) =A, 



(21) 



(22) 



where A is the Lagrange-multiplier used to introduce the condition (|10|) in the 
action S. So the first term of the right-hand side of (pi]) becomes 



SF 



5p{z) 



[Ps{z) - p w (z)]dz = 

= o, 



p s {z)dz - / p w {z)dz 



(23) 



in which we have used the fact that p w and p s , both satisfy (|10|) — perhaps with 
different integration regions, but it is not important, as one may define p to be 
zero outside the integration region and then extend the integration region to 
the whole real line. Next, consider the second term of (pl|). One can divide the 
integration region into three parts: 



i?2 

R 3 



z,Z G L s , 
z, Z G L w , 

z G L s and z' G L w or visa versa, 
where L s is the region in which p w > 1, and L w is the remaining region. We 
have 

Pw(z) - p s (z) ~ a, z G L s . (24) 
From this, it is seen that 



[pw(z) - p s (z)]dz - a 3/2 , 



(25) 



as the width of the region L s is of the order a 1 / 2 ©. Using © for p s and 
p w , it is seen that 

[p s (z)- Pw (z)]dz~a 3/2 . (26) 
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The width of L w is of the order a . This shows that on the average, 

Ps{z) - p w {z) ~ a 3/2 , z€L w . (27) 

Using these, it is seen that in Ri, the integrand in the second term of (|T]) is of 
the order a 2 , while the area of the integration region is of the order a. So the 
integral in this region is of the order a 3 . In i?2, the integrand is of the order 
a x a 3 / 2 , and the area of R2 is of the order a 1 / 2 . So the integral in this region 
too is of the order a 3 . Finally, in R3, the integrand is of the order a 3 / 2 x a 3 / 2 , 
while th area of R3 is of the order a . So, the integral in this region is of the 
order a 3 as well. One therefor deduces 

F s -F w ~a 3 . (28) 

Here we have used p'^(zo) 7^ 0. 

To find the critical behavior of the system, we must obtain F s — F w as a 
function of A — A c . To do so, consider a as a function of A and expand it 
around A = A c : 

a(A) = a(A c ) + a'(A c )(A - A c ) + ■ ■ ■ . (29) 
But a(A c ) = 0. So if a'(A c ) ^ 0, we have 

a~(A-A c ), (30) 

and from this 

F s - F w ~ (A - A c ) s , (31) 

that is, the model has a third-order phase transition. What we have proved, is 
that if Pw{zq) ^ and a'(A c ) ^ 0, then the order of the phase transition is 3. In 
fact, to prove this we have made an implicit assumption that ( p7j ) holds not only 
on the average, but at almost every point in L w . There might be very special 
cases where the difference between p s and p w is larger than a 3 but changes sign 
exactly in a such a manner that the integral of this difference is of the order 
a 3 / 2 . Also note that the whole reasoning is independent of the number of points 
at them p w attains its absolute maximum. 

Continuing to make this assumption, let's see what happens if either p'^(zo) 
or a'(A c ) is zero. As zq is the point which maximizes p w , the first nonzero 
derivative of p w at this point should be an even derivative. Let it be the (2m)- 
th derivative. Also, a should be increasing with respect to A at A c . This shows 
that its first nonzero derivative should be on odd derivative. Let it be the 
(21 — l)-th derivative. It follows that the width of L s is of the order a 1 ^ 2 ™), 
and from this 

F S ~F W ~ a 2+(1/m) . (32) 

One also has 

a~{A-A c ) 2l ~\ (33) 
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Combining these two, one arrives at 

F a -F w ~ (A-AJW-W+W^, 
where m and I are positive integers. 



(34) 



3 G(<p) = (j) 2k gYM 2 's 



For G{4>) = <p 2k (k an integer greater than 1), it has been shown that the density 
function p w (z), has a minimum at z = and two symmetric maxima at 

points z = ±akJjE, [(Tgl - is the parameter appearing in (|l6|), and can be 
evaluated from (fL7|) : 

r 2 fc (fc-l)! 



flfc _A(2fc-l)!! 
j/fc is independent of a& and is determined from 

k ~~ (2n-l)!! _ (n+1) 



0. 



Now it is possible to find the ^4-dependence of p(zo). Using (16) 

■ " 1 (2n- 1)!! 2n 2fe-2n-2 



Ptu(-z) — y <i; 

7T 



2™n! 



-a,, z 



(35) 



(36) 



(37) 



From this, 



Pw{zo) 



A2k 



2 k {k-l)\ 



fc-i 

£ 

n=0 



(2n-l)H 



2"n! 



-Vk 



{2k- 1)!! 

=: A&/(fc), (38) 
where /(fc) is independent of A. A c is determined by Pk(zo) = 1, which yields 

1 



A r = 



From this, 



«{£) = ( ^ 



[/(*)] 



.4 



2k ' 



1, 



(39) 



(40) 



which shows that the derivative of a with respect to A never vanishes. It is also 
easy to check that the p'^(z ) doesn't vanish either. Using p' w (z ) — 0, it is seen 
that 



Pw( z o) 



kAa 2k - s (2k - l)u 

7T \/l - U 2 



fe-2 



(2n-l)!! n k _ 



-(2fc-i) M 2fc - 3 + y v :"~ x; -\ 



2n-5 
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kAa 2k - 3 2k -1 ^ (2n-l)\\ 2fc _ 2 „_ 4 

v n— 

which is clearly negative. Here 

u:= Z A (42) 
o 

and © has been used. We conclude all gYM 2 's with G(0) = 2fc gYM 2 's 
exhibit third-order phase transition. 



4 The model G(<p) = (j) 2 + # 4 

We now consider G(<^>) — <p 2 + g(f) A , with g > so that in this G has a minimum. 
Using (|l6|), we have 

p w (z) = -Va 2 - z 2 (l + ga 2 + 2gz 2 ), (43) 

and from (p"7|), 

iA fl 2 + | 9 Aa 4 = l. (44) 

The shape of ^(z) depends on A. For 3ga 2 < 1, or A > 4g using (p4|), /^(z) 
has only one maximum at z = 0. For A < Ag, it has a minimum at z = and 
two maxima at ±z . We consider three distinct cases. 



4.1 g<A c /A 

In this case, p^, has two maxima for small areas. However, if the area exceeds 
4g < A c , p w will have only one maximum. As this happens before the critical 
area, the phase transition properties are governed by a one-maximum density 
function. The maximum of this function occurs at z = 0. Evaluating a from 
([fil), we arrive at 



(45) 




and from this 



a W = (3gMc) - 1 + ^ TTT ^ > 0. (46) 

V 27 ^ y/-l + y/l + 12g/A c 

It is also seen that 

^(0) = — (3 5 a 2 -l)<0. (47) 
na 

So the order of the transition is 3. 



4.2 g > A c /4 

In this case, for areas up to 4g > A c , the density function has two maxima at 
±z : 

Z ° = \fl (- 2+ ^l)' 
From this, one obtains p w {zo), and then 

It is also not difficult to check that p'w(zo) is negative. So in this case too, the 
order of the phase-transition is 3. 



4.3 g = A c /4 

In this case, for A < Ag — A c the function p w (z) has two maxima at ±zo and one 
minimum at z = 0. As A is increased, the difference p w {zo) — p w (0) becomes 
smaller, and at A — A c it becomes zero. For A > A c , p w (z) has only one 
maximum at z = 0. At A = A Cl the maximum of the density function becomes 
1 (and occurs at z = 0). Using this, the form of p w from ([f3"|), and ([ji]), it is 
seen that 

27 

Ac = ^\ (50) 
50 = T = 256 7r ' (51) 

and 

a c = g. (52) 

The interesting point for g — g is that p'^(0)\a=a c =0. In fact, one can see 
that 

„, , 3goAa(a + a e ) 

= -^(A-A c ) + -.., (53) 
while pL 4) (0) is negative. So for z 1; satisfying p w (zi) = 1, or 

7^(0)^ 2 + ^P { w 4) (0)zt = 1 - Pto (0), (54) 

one obtains 

N-a 1 ^ (55) 
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One also has 



from which it is seen that 



28 
27^ 



(A - A c ) + ■ ■ ■ , 



(56) 



|*l|~(A-4:) 1/4 . 



(57) 



This is an example of the case m = 2 and 1 = 1, discussed at the end of section 
2. It follows then that 

F s -F W ~{A-A c fi 2 . (58) 

The model G(<f>) = <fi 2 + (277r 2 /64)0 4 exhibits thus a phase transition of order 
5/2. 

5 Explicit evaluation of the order of the phase 



transition of the model G((f>) 



+ 



27tH 
64 



We now try to study the phase transition of the model G{cf>) — (f> 2 + go4> 4 
explicitly (g = 277r 2 /64). Using JTe| ) and (0), one finds that in the weak 
region 

FUA) = + ^g a 6 A + ^9 2 a s A, (59) 

. In the strong region, one uses the following 



where a can be found from 
ansatz for p s (z) 

Ps{z) = 



1, z€[-b,b] 

p s (z), z e [-a, -b] \J[b, a] 



(60) 



(as p w (z) has only one maximum for A > A c ). Using the standard method of 
solving the integral equations obtained (see [fl8| for more details) one arrives at 



F' S (A) = 



-M 



-N 



dX 



-M 



b \2 



X dX 



-b U(X) J_ b U(X) 



b a 



X 4 dX 



-go 



4 16 J J_ b U(X) 

+ -M 2 N / — — 

^8 2 )J_ h U{\) _ 

i r b x 4 dx _ r b x 6 dx 

+ 2 J_ b U{X)~ J_ b U(X)_ 
-2g A(~MN-±M 3 

1 N ^ + L M 2 N-—M £ 
8 16 128 



(61) 
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in which 



U(X) 
M 
N 



= V(a 2 -A 2 )(62-A 2 ), 
= a 2 + b 2 , 

= a 2 b 2 . (62) 



The parameters a and b are determined from following relations 

d\ 



b J\ 

A + go AM I — - - 0, (63) 



and 



b X 2 d\ 



2 MA + *> A [* M - N )-J_ b UiX) =1 - (64) 

As we are interested in phase transition, it is sufhcient to look at the solutions 
for A slightly greater than A c , which we obtain perturbatively. 
To do so, we use the following change of variables 

a =: a c (l — h), (65) 

and 

(b/a) 2 =: v, (66) 



in which a c is given by ( p2[ ) . Eliminating A from (|63|) and (^4|) , and solving the 
result for h, one obtains 

27 2 131 3 

h = v 2 H v 3 H , (67) 

224 6272 ' V ; 

from this, and (|63|), one obtains A in terms of v: 
A-A c 9 9 109 o 

— ; — - = — v 2 H v 3 H . 68) 

A c 28 784 v ' 

Note that A — A c is proportional to b 4 (in the leading order), as expected. 
Substituting these in (|6l|), one obtains F' a : 

F' s = 2 (36848 - 10836?; 2 + 6289t; 3 + ••■)• (69) 

From (fj~4|), one obtains a in terms of A. Substituting this in ([59]), and using 
( |68| ) as a definition for w, one obtains 

^ = 35721^ ( 36848 ~ 10836w2 " 4687w3 + •••)' ( 7 °) 

from which 

F' s {A)~F' w {A) = ^v 3 + ---. (71) 
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Using ([38|), one can write this as 

which proves that the order of the phase transition is 5/2, in agreement with 
©■ 
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